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ABSTRACT 

ry^ ■ Aims. We aim to investigate the long-term temporal evolution of the magnetic field of the solar-type star £ Bootis A, both from direct 

magnetic field measurements and from the simultaneous estimate of indirect activity indicators. 
i-C 1 Methods. We obtained seven epochs of high-resolution, circularly-polarized spectra from the NARVAL spectropolarimeter between 

2007 and 2011, for a total of 76 spectra. Using approximately 6,100 photospheric spectral lines covering the visible domain, we 
employed a cross-correlation procedure to compute a mean polarized line profile from each spectrum. The large-scale photospheric 
magnetic field of the star was then modelled by means of Zeeman-Doppler Imaging, allowing us to follow the year-to-year evolution 
of the reconstructed magnetic topology. Simultaneously, we monitored the width of several magnetically sensitive spectral lines, the 
radial velocity, the line asymmetry of intensity line profiles, and the chromospheric emission in the cores of the Ca II H and Ha lines. 
Results. During the highest observed activity states, in 2007 and 2011, the large-scale field of £ Bootis A is almost completely 
axisymmetric and is dominated by its toroidal component. The toroidal component persists with a constant polarity, containing a 
■ significant fraction of the magnetic energy of the large-scale surface field through all observing epochs. The magnetic topologies 

^ ' reconstructed for these activity maxima are very similar, suggesting a form of short cyclicity in the large-scale field distribution. 

, The mean unsigned large-scale magnetic flux derived from the magnetic maps varies by a factor of about 2 between the lowest and 

\^ . highest observed magnetic states. The chromospheric flux is less affected and varies by a factor of 1.2. Correlated temporal evolution, 

due to both rotational modulation and seasonal variability, is observed between the Ca II emission, the Ha emission and the width of 
l/^ 1 magnetically sensitive lines. The rotational dependence of polarimetric magnetic measurements displays a weak correlation with other 

activity proxies, presumably due to the different spatial scales and centre-to-limb darkening associated with polarimetric signatures, 
as compared to non-polarized activity indicators. Better agreement is observed on the longer term. When measurable, the differential 
rotation reveals a strong latitudinal shear in excess of 0.2 rad d~' . 
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^| 1. Introduction dwarfs (e.g. iBaliunas et alj Il995|: ISaar & Brandenburg 



... ... . TJM 

b. . lLockwood et al.ll2007t lOlah et alJl2009t iMetcalfe et al.ll2010l) . 



It is generally accepted that the regular succession of mag- More recent i Vj spectropolarimetric observations of cool dwarfs 

netic minima and maxima observed in the Sun and in many have become sufficiently accurate to enable t he direct detectio n 

cool stars is the result of an astrophysical dynamo, triggered of ma g ne tic fields on low-activity stars (e.g. iPetit et al1 l2008). 

by the combined presence of an outer convection layer and further expanding our insight into magnetic variability through 

stellar rotation (e.g. |Brandenburg & Subramanian||2005|). Using the ability to mon itor the long-term evolution of magnetic vec- 

chromospheric emission as a magnetic proxy, observations of torial topologies, instead of the disc-averaged chromospheric or 

stellar cyclicity have been conducted for several decades, re- photospheric fluxes 
vealing a richness in the temporal behaviour of solar-type 

For most targets, the time-base of spectropolarimetric obser- 



Send offprint requests to: A. Morgenthaler vatlons 18 stl11 restricted to a few years. This limited time span 

* Based on observations obtained at the Bernard Lyot Telescope can onl y offer a fragmentary view of magnetic cycles, except 

(TBL, Pic du Midi, France) of the Midi-Pyrenees Observatory, which for cycle periods much shorter than solar. A few dwarfs have 

is operated by the Institut National des Sciences de l'Univers of the already bee n observed to un dergo at least one global polarity 

Centre National de la Recherche Scientifique of France. switch (e.g. iPetit et alJ l2~009h . or to complete a full magnetic 
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cycle dFares et al ] 120091: iMorgentha leretalJ 1201 lh . Following 
the temporal evolution of large-scale stellar magnetic topolo- 
gies provides important constraints for numerical simulations 
of stellar dynamos, particularly now that the magnetic cycles 
of cool stars can be investigated through 3-D MHD simulations 
dBrown etalJl201d l201lt iGhizaru et alJl20Toh . 

Systematic comparisons between the temporal evolution of 
the large-scale magnetic field and other activity proxies (pho- 
tospheric or chromospheric) are still largely unexplored. Such 
studies are of interest because different measurable quantities 
related to magnetic activity can carry complementary informa- 
tion about the magnetic field generation in cool stellar objects, 
in particular through the different spatial scales to which they 
relate. As a step in this direction, we concentrate here on the 
solar-type star £ Bootis A. This active star is a main-sequence 
dwarf with an effective su rface temperature of 56 00 K and a 
surface gravity log g-A.65 (IValenti & F ischer 2001). As a visual 
binary system, the masses of £ Bootis A and B have been de- 
termined from astrometry to be 0.85 and 0.72 M respectively 
(|Wielenlll962h. The v ery high magnetic activity of £ Bootis A 
(Baliunas et al. 1995) has allowed early magnetic field detec- 
tions (Robinson et al. 1980) and is linked to a fast rotation period 
of 6.43 d (|Toner & Gravlll988l) . Given a low v sin i o f 3 kms" 1 
(iGrayl 11984 and a stellar radius of about 0.8 R Q dPetit et alJ 
2005), the short rotation period implies a low stellar inclination 
angle of about 28°. 

In this article, we simultaneously investigate the seasonal 
evolution of various activity proxies of £ Bootis A using (a) 
the large-scale surface magnetic topology, along with its short- 
term distortion through latitudinal differential rotation, (b) the 
Zeeman broadening of high Lande factor spectral lines, (c) the 
Ca II H and Ha core emission and (d) the radial velocity of in- 
tensity profiles, together with their asymmetry. We first describe 
the instrumental setup and spectropolarimetric time-series used 
in this study and the procedure employed to extract Zeeman sig- 
natures. We then detail the reconstruction of the large-scale mag- 
netic topology of the star at seven different epochs, followed by 
the extraction, from the same data sets, of a number of classic 
activity tracers. We finally discuss the results derived from our 
measurements. 



2. Instrumental setup, data reduction, and multi-line 
extraction of Zeeman signatures 

The data sets were collected with the NARVAL spectropolarime- 
ter at the Telescope Bernard Lyot (Observatoire du Pic du Midi, 
France) w hose instrumental setup is identical to the one de- 
scribed by iPetit et ail d2008l) . The echelle spectrograph has a 
resolution of 65,000 and covers the whole optical wavelength 
domain from near-ultraviolet (370 nm) to near-infrared (1,000 
nm) in one exposure, with 40 orders on the CCD. NARVAL pro- 
vides simultaneous recordings of the intensity spectrum (Stokes 
I) and one polarized spectrum (linear or circular). Here we deal 
with circular polarization (Stokes V parameter). Seven sets of 
spectra were obtained between 2007 and 201 1 (see Tableland 
|4]). They each contain between 7 and 19 spectra (in 201 1.07 and 
2008.09, respectively), obtained over 15 to 62 consecutive days 
(in 2007.59 and 2010.04, respectively). 

The raw data were processed using LibreEsprit, an automatic 
reduction softwa re developed for N ARVAL, based on the algo- 
rithm detailed bv lDonati et al.1 (Il997l) . For the reduced spectra of 
£ Bootis A, the maximum signal-to-noise ratio (S/N) is achieved 
at a wavelength of about 730 nm. The S/N value depends on the 



adopted exposure time and weather conditions above the Pic du 
Midi Observatory, with typical values varying between 450 and 
700. The S/N drops in the blue part of the spectrum, with a usual 
value of about 100 around the Ca II H&K lines. The situation is 
better in the near-infrared, with a S/N close to 500 in the neigh- 
bourhood of Fel@846.8404 used later in this study. 

The S/N of the polarized spectra is not sufficient to detect 
Zeeman signatures in individual spectral lines, even for photo- 
spheric lines with high Lande factors. To achieve a sufficiently 
low noise level, we used the reduced spectrum to calculate a 
single, cross-correlated photospheric line profile usin g the Least 
Square Deconvolution (LSD) multi-line technique (Donati et al. 
1997; Kochukhov et al. 2010). We employed a line-list of about 
6,100 spectral lines matching a stellar photospheric model for 
the spectral type of ^ Bootis A (G8V). This line -list is the same 
as the one previously used bv lPetit et alJ (12005). Thanks to this 
cross-correlation approach, the noise level of the mean Stokes V 
line profiles is reduced by a factor of about 30 with respect to 
the initial spectrum, so that the resulting noise level lies in the 
range 3.0 x 10~ 5 - 1.6 x 10 _4 / c , where I c denotes the continuum 
intensity. 

Fig.[T]shows the seven sequences of Stokes V LSD line pro- 
files corresponding to the seven observing runs. The rotational 
phases, indicated in the right-hand part of the plot and listed 
in Tab. [3] an d HI were calculated according to the ephemeris 
of Petit et al. (2005), who adopted a rotational period of 6.43 d 
dToner & Gravf l988) and a Julian date for zero rotational phase 
equal to 2,452,817.41. This rotation period is longer than the 
equatorial period derived from the modelling of the surface dif- 
ferential rotation (see Sect. 13.21 ). because it better represents the 
rotation period of the higher stellar latitudes that are expected to 
dominate the disc-integrated activity tracers, owing to the star's 
low inclination angle. The rotational modulation of most activ- 
ity proxies investigated in this study is more pronounced when 
using the adopted period of 6.43 d. 

Polarized features in the profile core, interpreted as Zeeman 
signatures, are detected for a large fraction of the observations, 
with amplitudes comfortably exceeding the noise level. It can al- 
ready be seen in Fig. [T] that the signal amplitudes are higher in 
2007.59 than in the following years. If we take a closer look at 
the shape of the signatures, we note that the majority of them 
are antisymmetric about the line centre (e.g. the rotational phase 
1.959 in 2007.59). A few of them are almost symmetric (phase 
3.203 in 2007.59, 4.04 in 2008.09, and 10.172 in 2010.04). We 
note also that the sign of the signatures (symmetric or antisym- 
metric) is the same throughout the observations. 

3. Photospheric magnetic field 

3. 1 . Longitudinal magnetic field 

Using the LSD line profiles in the Stokes parameters / and V, we 
can derive the longitudinal component B/ of the surface mag- 
netic field, as averaged over the visible hemisphere of the star. 
To do so, we used t he centre-of-gravity method described by 
iRees & Semell (Il979l) . through the equation : 

,, JvV(v)dv 

B,(G) = -2.14 x 10 11 J — , (1) 

Aogc J(l-/(v))</v 

where v (km.s~ l ) is the radial velocity in the stellar rest-frame, 
/lo (nm) the mean wavelength of the line mask used to compute 
the LSD profile (538 nm in our case), g the mean effective Lande 
factor of the line list (equal to 1.21), and c (km.s^ 1 ) the velocity 
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Fig.l. Normalized Stokes V profiles (V/Q of £ Boo A for 2007.59, 2008.09, 2009.46, 2010.04, 2010.48, 2010.59 and 2011.07 
(from left to right and top to bottom), after correction for the mean radial velocity of the star (the radial velocity values are listed in 
Tab.|2]i. Black lines represent the data and red lines correspond to synthetic profiles of our magnetic model. Successive profiles are 
shifted vertically for display clarity. Rotational phases of observations are indicated in the right-hand part of the plot and error bars 
are illustrated on the left side of each profile. 



of light in vacuum. The integration window covers a velocity The B[ values are listed in Table [5] and [6] and the aver- 
range of +16 kms~' around the line centre. aged values for each epoch are recorded in Table Q] We note 

that for Stokes V profiles with the same amplitude, longitudinal 
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Table 1. Magnetic quantities derived from the set of magnetic maps. 



frac. year 


Timespan 


n0 


" hong 


"mean 


pol. en. 


dipole 


quad. 


oct. 


axi. 




dQ. 


xi 


(2000+) 


(d) 




(G) 


(G) 


(% tot) 


(% pol) 


(% pol) 


(% pol) 


(% tot) 


(radd- 1 ) 


(radd" 1 ) 




07.5872 


14.99 


10 


8.7 + 7.1 


69 + 27 


17 + 3 


71+2 


13+ 1 


9+1 


83 + 3 






4.8 


08.0881 


28.04 


19 


4.6 + 3.1 


30 + 8 


58 + 6 


41+7 


15+ 1 


15 + 2 


56+1 


1.13 + 0.01 


0.38 ± 0.02 


1.78 


09.4572 


38.02 


13 


8.3 ± 6.4 


47+11 


35 + 9 


41+7 


20+1 


19 + 3 


69 + 3 


1.27 + 0.01 


0.57 ± 0.03 


1.37 


10.0403 


61.85 


9 


4.1+5.3 


38 + 9 


32 + 9 


29 + 20 


9+1 


8 + 1 


29 + 6 






1.3 


10.4795 


37.98 


10 


8.6 + 3.9 


36+12 


62 + 2 


50 + 6 


13 + 3 


7 + 2 


43 + 8 


1.055 + 0.005 


0.67 ± 0.01 


2.5 


10.5945 


23.93 


9 


11.2 + 4.2 


46+17 


13 + 6 


48 ± 14 


18 + 4 


13+4 


95 ± 1 


1.09 + 0.03 


0.27 ± 0.05 


1.64 


11.0657 


17.99 


7 


8.4 + 3.1 


43 + 20 


18 + 4 


77 + 3 


14+1 


4 + 3 


85 + 2 






2.19 



Notes: As a function of the mean fractional year of each run, we list the timespan of successive data sets, the number of available rotation phases, 
the average longitudinal field (with its associated dispersion), the mean unsigned magnetic field strength (S racan ), the fraction of the large-scale 
magnetic energy reconstructed in the poloidal field component, the fraction of the poloidal magnetic energy in the dipolar (I = 1), quadrupolar 
{E = 2), and octopolar {t = 3) components, and the fraction of energy stored in the axisymmetric component (m = 0). We then list the differential 
rotation paramet ers fi cc| and dQ and the^ of the magnetic models. Error bars on the ZDI magnetic quantities are derived using the method of 
iPetit et all d2008h . 



Table 2. Mean values of the activity tracers. 



fractional year 


v r 


Mean velocity spans 


Mean line widths 


NcallH 


log R' HK 




(2000+) 


(km s~') 


(kms" 1 ) 


(krns" 1 ) 










07.5872 


1.86 + 0.02 


0.03 ± 0.02 


15.09 + 0.28 


0.4434 ± 0.0080 


-4.34 + 0.01 


0.3588 ± 


0.0020 


08.0881 


1.74 + 0.05 


0.04 ± 0.02 


14.59 ± 0.28 


0.4198 + 0.0079 


-4.38 ± 0.01 


0.3552 ± 


0.0015 


09.4572 


1.88 + 0.03 


0.03 + 0.01 


14.34 ± 0.23 


0.4207 ± 0.0047 


-4.38 ± 0.01 


0.3532 ± 


0.0012 


10.0403 


1.78 + 0.04 


0.03 ± 0.01 


13.79 + 0.20 


0.4004 ± 0.0060 


-4.40 ± 0.01 


0.3513 + 


0.0010 


10.4795 


1.93 + 0.03 


0.04 ± 0.02 


14.15 + 0.16 


0.4032 ± 0.0059 


-4.40 ± 0.01 


0.3512 + 


0.0011 


10.5945 


1.94 + 0.02 


0.05 + 0.01 


13.96 ± 0.28 


0.4027 + 0.0114 


-4.40 ± 0.02 


0.3503 ± 


0.0016 


11.0657 


1.83 + 0.02 


0.03 + 0.01 


14.52 + 0.33 


0.4281 ±0.0108 


-4.37 + 0.01 


0.3573 ± 


0.0025 



Notes: we list the radial velocity of the star, the bisector span of Stokes I mean profiles, the width of the Fel@846.8404 line, the Ca II emission 
index (along with its \og(R' HK ) counterpart), and the Ha emission index. All average quantities are given with their dispersion. 



fields derived using the centre-of-gravity technique are a maxi- 
mum when the Stokes V signature is anti-symmetric about the 
line centre. This is often not the case in our observations, be- 
cause £ Bootis A has a large-scale toroidal field component (see 
lPetitetal.ll2005l and Sect. EO l. If a profile has a higher red- 
shift for the zero-crossing velocity, for example, phase 6.148 in 
2010.48, this can result in a lower Bi value, in spite of a high 
Stokes V amplitude. 

The rotational dependence of B\ (Fig. [3] and HJ is visible 
at several epochs, and is especially obvious in 2007.59 and 
2010.59. The phase dependence is generally more complex than 
a sine curve, indicating that the surface magnetic geometry is not 
limited to a simple dipole. A long-term trend is also observed in 
the averaged Bi values (Fig. [6), with a strength of 8.7 G (with a 
dispersion of 7.1) in 2007.59 and a much lower average value of 
4.6 G (with a dispersion of 3. 1) in 2008.09. This global decrease 
can readily be seen in the decreased amplitude of the Stokes V 
profiles between the two epochs. Most of the later measurements 
have a similar or higher B\ value to that in 2007.59, with the ex- 
ception of 2010.04, which has a value similar to that in 2008.09. 

3.2. Magnetic maps 

The rotational modulation observed in the line-of-sight projec- 
tion of the magnetic field (Fig. [TJl indicates a complex surface 
distribution of magnetic fields. We can model this complex field 
topology using a tomographic approach. 

We employed Zeeman-Doppler Imaging (ZDI) to recon- 
struct the surface distribution of the magnetic vector from the 
time-series of Stokes V LSD profiles at eac h observing epoch. 
This method, first propos ed by Semel ( 1989), was implemented 
and tested by IDonati & Brownl d 19971) . The version of the ZDI 
code used here assumes that the field geometry is projected onto 



a spherical harmonics frame (IDonati et al.ll2006h . In this inver- 
sion procedure, the time series of Stokes V LSD profiles is com- 
pared to a set of synthetic Stokes V line profiles computed for the 
same rotational phases as the observed profiles. Synthetic Stokes 
profiles are calculated from an artificial star whose surface is di- 
vided into a grid of pixels. Each surface pixel is associated with 
a local Stokes I and V profile. Assuming a given magnetic field 
strength and orientation for each pixel, local Stokes V profiles 
are calculated under the weak-field assumption, where Stokes V 
is proportional to g.A^.B\\.dl /dA, with Aq representing the aver- 
age wavelength of the LSD profile, B\\ the line-of-sight projec- 
tion of the local magnetic field vector and g the effective Lande 
factor of the LSD profile. 

We furthermore assumed that there are no large-scale bright- 
ness inhomogeneities over the stellar surface, so that all syn- 
thetic Stokes I profiles are locally the same everywhere. The lim- 
itation of this assumption will be illustrated by the observed vari- 
ations in radial velocity and profile bisectors (Sect. l4~2l i. provid- 
ing us with evidence that Stokes I profiles do vary as a function 
of the rotational phase, due to starspots or plages. The restric- 
tion of the reconstructed magnetic topology to the global-scale 
component of the surface magnetic field is probably limiting the 
consequences of these model inaccuracies, which should affect 
mainly smaller, unresolved spatial scales. In the present model, 
the Stokes I profiles are assumed to possess a Gaussian shape, 
with a depth and width adjusted to achieve a good match be- 
tween synthetic and observed Stokes I line profiles. We assumed 
a projected rotational velocity (v sin z) of 3 km s~' and an incli- 
nation angle equal to 28°, both of these values were previously 
used in th e for ward modelling of the field topology made by 
IPetit et ail (120051) . We also assumed the limb-darkening to be 
linear with fi = cos(fl), where 9 is the limb angle, with a coeffi- 
cient equal to 0.75. We finally restricted the spherical harmonics 
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Fig. 2. Magnetic maps of £ Bootis A, derived from 2007.59, 2008.09, 2009.46, 2010.04, 2010.48, 2010.59, and 201 1 .07 observations 
(from left to right and top to bottom). For each data set, the three charts illustrate the field projection onto one axis of the spherical 
coordinate frame with, from top to bottom, the radial, azimuthal, and meridional field components. The magnetic field strength is 
expressed in Gauss. 



expansion to £ m . dx < 10, since no improvement to the fits, be- 
tween modelled and observed LSD profiles, is noticed for values 
of (max greater than 5. 

Because each data set was collected over several weeks, we 
assumed that the magnetic geometry might be distorted by lati- 
tudinal differential rotation over the timespan of the data collec- 
tion. We therefore included a two-parameter differential rotation 
law in our stellar model, with the form 

Q(/) = C2 eq - dQ. sin 2 I , (2) 

where Q(7) is the rotation rate at latitude /, Q eq the rotation rate 
of the equator and c/Q the difference in rotat ion rate between the 
poles and equator. Following the method of [Petit et alj d2002l) . a 
grid of magnetic inversions was calculated for a range of values 



of the quantities £2 eq and dQ.. The values listed in Tab. Q] corre- 
spond to the x 2 minimum in the parameter plane, whenever this 
minimum exists and is unique in the scanned Q. eq -dQ. area. 

One first limitation of this indirect imaging procedure is the 
roughness of the underlying stellar model, which, combined with 
the sparse phase sampling and uneven S/N ratio, may be the 
source of inaccuracies in the reco nstructed magnetic geometry 
(see, e.g.. iDonati & Brownl [l997). One other limitation of the 
maximum-entropy algorithm is the absence of error bars on the 
resulting maps. To limit the consequences of these two restric- 
tions as much as possible, we do not discuss here the finest de- 
tails of the magnetic topology (i.e., individual magnetic spots), 
but rather concentrate on a set of quantities derived from the 
largest spatial scales of the field geometry (e.g., its low-order 
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Fig. 3. Four sets of six subplots, each set corresponding to one epoch (from top to bottom, 2007.59, 2008.09, 2009.46, and 2010.04). 
For each set, the left subplots contain, from top to bottom, the A^-index (filled circles and left axis) and NcaiiH -index (crosses and 
right axis of the top panel of each set), the Fel@846.8404 magnetic line widths (km.s~ l ), B/ (Gauss). The right subplots, from top 
to bottom, correspond to the velocity spans (km.s~ l ) obtained from the profile bisectors (Sect 4.2), the radial velocity (km.s~ l ), and 
the mean unsigned magnetic strength of the total (full line), radial (dashes), azimuthal (dot-dashes) and meridional (dots) magnetic 
field components extract from the ZDI maps. Error bars are not included whenever they are smaller than the symbol size. 



multipolar expansion), as listed in Tab. Q] As an attempt to eval- 
uate the uncertainties on these global magnetic quantities, we 
then reproduce the approach of iPetit et alJ (120081) and compute 
several magnetic maps, each of which is calculated using dif- 



ferent values of the input parameters of ZDI (within the error 
bars on individual parameters). The dispersion of the resulting 
magnetic values are c onsidered as error bars. As for other solar- 
like stars presented bv lPetit et alJ d2008l) . we generally note that 
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Fig. 4. Same as Figure for 2010.48, 2010.59 and 2011.07. 
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the resulting dispersion is dominated by the uncertainty on the 
inclination angle of £ Bootis A. 

The fits from our tomographic modelling are illustrated in 
Fig. [T] and the magnetic maps are given in Fig. [2] The mod- 
elled data were fitted to the observed ones with a reduced y 1 
(x\ hereafter) of between 1.3 and 4.8 (see Tab. [TJ, showing that 
the fit accuracy never reached the noise level. The fit is not sig- 
nificantly improved by simple changes in the local line profile 
model, e.g. by the replacement of the Gaussian local Stokes I 
profile by a Lorentzian one (in an attempt to slightly improve the 
fitting of the profile wings), or the implementation of an ad-hoc 
asymmetry in loc al Stokes V profiles, as previously proposed by 
iPetit et ail (120051) . A possible cause for the r epeated mismatch 
may be the residual effect of blended lines dKochukhov et al.l 
2010) which could be larger than the noise level. Another likely 
cause is the intrinsic evolution of the magnetic topology during 
data collection, if this evolution is not entirely caused by to lat- 
itudinal differential rotation, or if the latitudinal shear is follow- 
ing a different law than the simple formula used here. Fast evo- 
lution of the magnetic geometry is observed between 2010.48 
and 2010.59, but the field configuration seems to be more stable 



at other epochs, for instance in 2010.04, where a 62 d timespan 
produced a magnetic inversion with^= 1.3. 

The limited lifetime of magnetic tracers is especially of con- 
cern in the measurement of surface differential rotation. In the 
ideal situation of a magnetic topology progressively distorted by 
a surface shear, the Xi lands cape in the Q ea - dQ. plane takes the 
shape of a 2-D paraboloid (IPetit et al]l2002l) . This is never the 
case here, except for 2008.09. At other epochs, iso-^ contours 
are sometimes not ellipsoidal (2009.46, 2010.48, 2010.59), de- 
limit at least two local ^ minima of similar depth (201 1.07), or 
display no minimum (2007.59 and 2010.04). Therefore the ap- 
parent temporal variations of the shear shown in Table Q] should 
be taken with caution, even when the observed differences are 
above the statistical error bars. In particular, we note that modi- 
fying the number of profiles used to model the surface shear (for 
instance, by removing the first or last profile of a time-series) of- 
ten has a significant impact on the measured dQ., suggesting that 
intrinsic magnetic variability is affecting some of our differential 
rotation measurements. Whenever a unique Xr minimum was not 
identified, the differential rotation parameters of 2008.09 were 
adopted in the inversion process. 
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The magnetic maps show the reconstructed geometry of the 
large-scale photospheric field. Over the years, the most recog- 
nizable feature is the azimuthal field component, organized in 
a fragmented or complete ring encircling the rotation axis at a 
latitude of about 40°, showing a large-scale toroidal magnetic 
component at the stellar photosphere. The toroidal component 
contains most of the magnetic energy, except in 2008.09 and 
2010.48 (Table [TJ, and its dominant polarity is constant for all 
data sets. The poloidal magnetic component is dominated by a 
dipole, with the dipole tending to host a larger fraction of the 
magnetic energy when the toroidal field component is strongest, 
i.e., in 2007.59 and 2011.07. The fraction of the axisymmetric 
field (fraction of magnetic energy stored in spherical harmon- 
ics modes with I = 0) is also evolving in correlation with the 
strength of the toroidal field. 

3.3. Zeeman broadening 

To investigate the photospheric magnetic field of £ Bootis Ain 
even greater, we looked at the broadening of several spectral 
lines with high magnetic sensitivity. T his approach, success - 
fully applied to £ Bootis A in the past (Robins on et al.| [l980), 
is complementary to the Zeeman-Doppler Imaging strategy, be- 
cause Zeeman broadening is sensitive to the total magnetic flux, 
while the polarized Zeeman signatures carry selective informa- 
tion about the large-scale component of the magnetic field (ow- 
ing to signal cancellation for close-by regions of opposite mag- 
netic polarity). Here, we monitored the temporal fluctuations of 
the Zeeman broadening. To achieve this goal, we adopted a very 
rough tracer of the magnetic broadening by simply following the 
time evolution of the line width. To do so, we first used a cubic 
spline algorithm to modify the wavelength sampling of the line 
to obtain a grid of spectral bins offering points of equal intensity 
on both wings of the line profile. We then evaluated the width at 
various depths in the line, and finally retained the depth (0.33/ c 
above the normalized intensity of the line center) at which the 
rotational modulation is the most obvious. The measured width 
is therefore not a standard FWHM. 

We have mainly focused on the Fel@846.8404 line, which 
offers a high Lande factor (g = 2.493) and an infrared wave- 
length, which both enhance the Zeeman broadening effect. Other 
advantages of using this specific line are the relatively high S/N 
obtained in this part of the NARVAL spectra (close to 500) and 
the clean surrounding continuum, which is mostly free from tel- 
luric or photospheric lines, which ensures accurate continuum 
normalization. The line width measurements are listed in Table 
[5] and plotted in Fig. [3] and [4] A complex phase dependence of 
the line width is obvious at several epochs. A correlated evolu- 
tion is sometimes observed with B\ (e.g. in 2010.48), but this 
occasional similarity cannot be taken as a general trend, as il- 
lustrated by the anti-correlation observed only about a month 
later, in 2010.59. In addition to the rotational modulation in line 
width, a long-term trend is observed, with a significant decrease 
in the line width between 2007.59 and the following data sets. 
This variability is larger than the rotational variability observed 
at individual epochs (Pig. [6j>. This decrease agrees with the trend 
in decreasing magnetic field strength seen in the magnetic maps 
(Tab.Q). 

The line width was also monitored using two other magneti- 
cally sensitive lines (Fel@549.75, g = 2.26 and Fel@550.68, g = 
2.00). The variations described above are consistently recovered 
using these lines, although the correlation with chromospheric 
emission (see Section 14.11 ) is generally worse than for the in- 
frared line. We also studied three lines with a weak Lande factor 
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Fig. 5. Ca II H line for August 2 2007 (solid line) and for January 
18 2010 (dashed line). 



(Fel@481.78, g = 0.48, Fel@592.74, g = 0.43 and Fel@840.14, 
g = 0.5 1 ) which, as expected, did not exhibit any significant tem- 
poral evolution in their width. 

4. Other activity tracers 

4. 1 . Chromospheric activity 

Using our Stokes I spectra we monitored the evolution of the 
chromospheric flux between each set of observations and also 
across the stellar rotation cycle. A simple comparison of the core 
of the Ca II H line for two different years (Fig. [5]) reveals the 
variability of the chromospheric activity, with the core emission 
stronger in 2007.59 than in 2010.04. We constructed two emis- 
sion indexes to quantify the emission changes, using the Ca II H 
and Ha lines, respectively. Our method is described below. 



4.1.1. Nr 



-index 



Before computing the index itself, the continuum normaliza- 
tion in the spectral region around Ca II H&K was adjusted. The 
very dense distribution of photospheric spectral lines in this part 
of the spectrum prevented the standard reduction pipeline of 
LibreEsprit from defining a reliable continuum level, resulting 
in local normalization inaccuracies of the order of 20%. To im- 
prove the situation, we took as a refere nce a synthetic norma lized 
spectrum from the POLLUX database (iPalacios et al.l2 010) with 
T e ff and log(g) values close to those of £ Boo A. First, we inter- 
polated the synthetic spectrum on the NARVAL wavelength grid, 
we then defined by hand a number of reference points around 
the Ca II H&K lines and used them to estimate the local ratio 
between the NARVAL spectrum and the synthetic spectrum. The 
series of ratios were then fitted by a fourth-orderpolynomial and, 
finally, the region of interest in the NARVAL spectrum was di- 
vided by the polynomial. We achieved a good continuum level 
using this method, without the need to rotationally broaden the 
synthetic spectrum. This confirms that rotational broadening can 
be ignored in the process, at least for a star like £ Bootis A that 
has a low v sin ;. 

Afterwards, w e calculated a Nr.ni/ H-index following the 
method outlined in iDuncan et alj (I1991I) . who defined S -values 
from Mount Wilson observations. We used two triangular band- 
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Fig. 6. Top three rows: long-term evolution of the average values (calculated over individual observing runs) and dispersion (vertical 
bars) of the activity proxies presented in Fig.[3]and[4] Bottom two rows: temporal evolution of magnetic quantities listed in Tab.Q] 



passes H and K with a FWHM of 0. 1 nm to determine the flux in 
the line cores. Two 2 nm-wide rectangular bandpasses R and V, 
centred on 400.107 and 390.107 nm, respectively, were used for 
the continuum flux in the red and blue sides of the H and K lines. 
Although results from this adapted S index are generally good, 
one concern was the location of Ca II K in an order overlap in the 
2007.59 data set (later on, the order overlap was slighlty shifted, 
so that the line core fell outside the overlap region). This unfor- 
tunate position in the spectrum generates additional difficulties 
in obtaining a reliable continuum normalization and increases 
the uncertainties in the computed values (since the S/N ratio de- 
creases rapidly towards the edge of an order). After multiple tri- 
als, we chose to use only the Ca II H line becau se it gave the 
most r eliable results. Finally, in the same way as Wright et al. 
(2004) defined an L-index for Lick observations, we defined a 
Ncaim-index for our NARVAL spectra as follows: 



_ H 

A 1 Cairn — — 
K 



(3) 



where H and R are the same as described above. 

To match the Mount Wilson values, we transformed the 
Ncaim-index as 



(H\ 2 H 



+ 7 



(4) 



where a, ft and y are relative weights to be determined. To esti- 
mate them, we chose 3 1 cool stars simultane ously present in the 
NARVAL archive and in Wright et al. (2004), and ensured that 
all selected NARVAL observations had an S IN above about 100 
around the Ca II H line. We then performed a least-squares fit- 
ting between the chromospheric activity values of the two stellar 
samples and found an optimal value of -0.972 for a, 1 .803 for 
/3 and -0.051 for y. Using the result of this calibration, some 
residual scatter is still observed between our measurements and 
Wright's values. These differences are probably due to the dif- 
ferent dates of observation between the two studies, bearing in 
mind that the chromospheric emission of cool active stars dis- 
play long-term fluctuations. The observed scatter, which is not 
associated to any systematic bias, is similar to that between 
Wright's estimates and older Mount Wilson measurements. 

The NcaiiH -index is shown in Fig. [3] and [4] and the mean 
values for each set appear in Table [2] Random errors are about 
10~ 3 for individual observations. Occasional repeated observa- 
tions of £ Boo A during a single night (twice each night on 2008 
Jan. 18 and Jan. 23) enabled us to obtain another estimate of 



uncertainties, assuming that the chromospheric activity is fairly 
constant over the few minutes that separate successive spectra. 
The typical difference in the NcaiiH -index between these close- 
by observations is about 4 x 10~ 3 . 

A rotational modulation in NcaiiH is visible in most data sets, 
and is particularly evident in 2008.09. A longer-term evolution 
is also seen (Fig. [6), with differences between the years being 
larger than the fluctuations observed during a single rotation cy- 
cle. In individual data sets, a correlation between the rotational 
modulation of NcaiiH and the width of Fel@846.84 is some- 
times observed (e.g., 2008.09) but is much less pronounced at 
other epochs (e.g., 2009.46). A possible reason for this partial 
mismatch is the different centre-to-limb behaviour of the two 
magnetic field tracers. Stokes I for Fel@846.84 will be affected 
by changes in the central, unsplit n component of the line (in- 
creasing towards the limb) as well as changes in the split <x 
components (decreasing towards the limb). The combination of 
the two types of Zeeman components gives a different centre-to- 
limb variation than for the chromospheric lines, for which a limb 
brightening is expec ted if there are fac ulae/plages contributing to 
the Ca II emission (Orti z et al.ll2002l) . The correlation between 
NcaiiH and Fel@846.84 is more easily seen on the longer term 
(Fig.0. 

4.1.2. JV# a -index 

In the same manner as for NcaiiH, we defined a quantity to mea- 
sure the variab ility in Ha. W e used the same rectangular band- 
passes that iGizis et al.l ([2002) defined around the Ha line. The 
red-side continuum is taken to be between 656.62 and 656.84 
nm, the blue-side continuum between 655.77 and 656.0 nm, and 
the line core between 656.10 and 656.46 nm. By defining the 
flux values in the continuum bandpasses as C re d and Cuue and 
the flux in the line as Fn a , our index was constructed as follows: 



N, 



'Ha 



Ha 



C re d + Cblu 



(5) 



Note that the higher S/N around Ha, as well as the better 
continuum normalization in this spectral region, allowed us to 
skip the renormalization procedure adopted for Ca II H. The re- 
sults obtained for each year are given in Figs. [3] and |4] and the 
mean values are listed in Table [2] Random errors are about 10~ 4 
for a single observation. Using the two nights for which repeated 
observations are available, as we did for NcaiiH, the empirical 
uncertainty was measured to be around 7 x 10~ 4 . 
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Fig. 7. Correlation between NcaiiH and Nh„ (top), and between 
the widths of the Fel@846.84 line and N C aim (bottom). 



Most conclusions drawn for the NcaiiH index are also valid 
for Nho, given the correlation of the two indices illustrated in 
Fig- 13 Again, rotational variability is observed for most epochs, 
as well as a year-to-year trend. The changes seen between the 
observing epochs are higher than the rotationally induced varia- 
tions. 



4.2. Radial velocities and profile bisectors 

Even for the low Doppler broadening of spectral lines suffered 
by £, Bootis A (vsin i= 3km s _1 ), surface inhomogeneities (cool 
spots or plages) can be expected to modify the shapes of the 
Stokes I LSD profiles. During its transit across the visible hemi- 
sphere of the star, a cool spot changes the depth of one wing 
of the profile when it gets closer to the stellar limb (so that its 
spectral counterpart gains significant Doppler shift). This effect 
may be detected as a change in the radial velocity of the line 
centroid, or as a modification of the profile bisect or. The typical 
profil e asymmetry produced by convective flows (IToner & Gravl 
1988) may also be locally modified by strong magnetic fields, 
resulting in a temporary change of the profile bisector. In this 
section, we propose to investigate these effects for £ Bootis A. 

Radial velocities (RV hereafter) were calculated from our 
Stokes I LSD profiles by applying a Gaussian fit to the line pro- 



files and using the centroid of the resulting Gaussian function as 
a radial velocity estimate. Using a (symmetric) Gaussian fit may 
be thought inadapted for asymmetric Stokes I LSD profiles ob- 
served for ^ Bootis A. We note, however, that the level of asym- 
metry is generally low compared to radial velocity fluctuations 
derived below (Table [2J. Detailed measurements can be found 
in Table [5] and phase-folded results are shown in Figs. [3] and [4] 
The typical error for a single measure of RV, using N ARVAL in 
polarimetric mode, was estimated by Mou tou et al.l (120071) and 
is of the order of 0.015-0.030 kms -1 . Several data sets display 
rotationally modulated variations in RV, with amplitudes some- 
times in excess of 0.1 kms -1 (see for instance 2009.46). Most 
of the time, this activity proxy is acting differently from other 
tracers. 

We also used the Stokes I LSD profiles to construct pro- 
file bisectors. After adjusting the sampling of the intensity line 
profile with a cubic spline p rocedure, we comput ed the bisec- 
tor using a method similar to IToner & Gravl d!988l) . The phase- 
dependence of the bisector shape is plotted in Figure [8] for the 
observations of 2008.09. We note that the bisectors never dis- 
play the typical redshift near the conti nuum that they exh ibited 
most of the time i n previous studies of lToner& Gravl d~1988) and 
iPetit et all d2005l) . This shape difference is due to the spectral 
range (covering the whole visible domain) used to derive the 
LSD profiles. Using various sub-sets of our line-list to compute 
LSD profil es in narrower spectral regions, we note that the shape 
reported bv lToner & Gravl ( fl9 88) is typical of the red- and near- 
infrared domains, while a mirror shape (with a blueshift close to 
the continuum) is observed for the bluest parts of the NARVAL 
spectral range. 

We note that strong and weak individual spectral lines 
(which are expected to have different bisector shapes) are mixed 
together in the computation of a LSD profile. If the resulting 
bisector shape is st i ll a tracer of photosp heric convective flows 
dGravlll980lll98U Il982t iDravind fl987h . any direct compari- 
son with individual line bisectors shou ld be taken with ca ution. 
When adopting the exact line list of IToner & Gravl dl988l) to 
compute LSD profiles (eight lines close to 625 nm), we indeed 
derived bisectors similar to theirs, but here again the comparison 
is of limited relevance, since it might be affected by the differ- 
ence in spectral resolution and time of observation. 

To quantify the changes in the bisector shape, we calcu- 
lat ed the velocity span a ssociated to each bisector, as defined 
by IToner & Gravl d 1 988b . Here, we take the difference in ra- 
dial velocity between a point near the top of the line profile 
(I/I c = 0.95) and a point near the bottom (the I/I c minimum 
value plus 0.01). Uncert ainties in the bisect or span were derived 
using the approach of iPovich et al.l (|2001), yielding error bars 
of 20-30 ms _1 . Repeated measurements within a few minutes, 
obtained at two occurences in 2008.09, provided us with vari- 
ations of the velocity span within the calculated uncertainties. 
The bisector fluctuations are shown in Figures [3] and |4l with nu- 
merical quantities listed in Tab. [5] Variations with the rotational 
phase are visible for most epochs of observation, and sometimes 
follow a complex pattern (e.g. 2008.09). Most of the time, no 
obvious correlation can be found with other magnetic or activ- 
ity proxies. In contrast to most other tracers, a long-term trend 
of the bisector span is not detected above its typical dispersion 
level (Fig. [6]). 

5. Discussion and conclusions 

We have used spectropolarimetric observations to derive a time- 
series of magnetic field maps and activity proxies covering al- 
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Fig. 8. Bisectors of Stokes I LSD profiles as a function of the rotational phase for 2008.09. Dashes represent observations taken a 
few minutes apart from the overplotted solid-line bisector. 



most four years in the life of £ Bootis A. In doing so we com- 
bined observational approaches that are often investigated sepa- 
rately. 



5. ) . Magnetic topology 

In spite of its rapid rotation rate, the spectral lines of £ Bootis A 
are not significantly affected by Doppler broadening. The low 
v sin i of the star is due to the combined effect of a low inclination 
angle and a relatively small radius. When applied to stars with 
low rotational broadening, ZDI selectively reconstructs only the 
largest spatial scales of the surface magnetic field. The recon- 
structed magnetic maps of £ Bootis A show a surface-averaged 
field strength of the order of 30-100 G. T hese values are si gnifi- 
cantly lower than the previous estimate of lPetit et alJ (120051) de- 
rived from observations taken during the summer of 2003, both 
for the poloidal and toroidal field components. This difference is 
possibly the consequence of a long-term magnetic trend, how- 
ever, we emphasize that the instrumental setup and modelling 
methods are di fferent in the two studies (direct modelling in 
iPetit et al.l2005l inverse tomography in the present study), which 
may account for part of the apparent decrease. Both methods 
agree on the existence of a strong toroidal field component on 
£ Booti s A, which is do minating the energy budget in the re- 
sults o f lPetitetal] (12005). and also in most of the ZDI magnetic 
field reconstructions in the present study, except in 2008.09 and 
2010.48 (for which the poloidal field component dominates the 
surface magnetic energy). 

The properties of the large-scale field of £ Bootis A can also 
be compared to that of other ZDI studies of solar-type dwarfs, as 
long as the selected stars also exhibit a low v sin i, so as to en- 
sure that all stars used for the comparison are affected by a same 
low-pass spatial filtering of their magnetic geometry. As far as 
mass is concerned , the closest available neighbour is HD 189733 
dFares et al.ll2010h . The average surface field strength recovered 
for this star is lower than that derived for £ Bootis A, which is 
likely a consequence of the slower rotation rate of HD 189733 
(with a rotation period of 12 d). In spite of this difference, both 
stars have a significant fraction of their surface magnetic energy 
reconstructed in a toroidal field component. 

A small number of main-sequence stars of higher mass 
can also serve as comparisons. The average field strength of 
£ Bootis A is consistent with the large-scale field strength mea- 
sured for 1 M H stars with rotation periods slower than about 
10 d (IPetit et al.ll2008l) . The prominent toroidal component of 



the large-scale magnetic geometry, observed for each year ex- 
cept for 2008.09 and 2010.48, agrees with the assertion, derived 
from both observations and numerical simulations of stellar dy- 
namos, that surface magnetic geometries become predominantly 
toroidal for solar-type dwarfs with rotation periods as short as a 
few days, whereas a ma inly poloidal surface field is a com mon 
feature of slow rotators (Petit et al. 2008; B rown et al 1 l2Ol0h . 



5.2. Magnetic variability 

If £ Bootis A seems to obey the above-mentioned trends, the 
temporal fluctuations of its magnetic properties demonstrate that 
magnetic variability will induce some scatter around any average 
behaviour, because the global field strength, the fraction of en- 
ergy in the toroidal field, and the field axisymmetry all vary with 
time. While relative fluctuations observed in the chromospheric 
emiss ion are limited to about 20% (although Baliuna set alJ 
1995, suggest that longer-term variations of up to 60% can be 
observed), simultaneous variations by a factor of about 2 are ob- 
served in the mean large-scale field strength, and reach a factor 
of about 5 in the fraction of energy reconstructed in the axisym- 
metric field component and a factor of 3 in the fraction of energy 
stored as a toroidal component. Assuming that these sharp vari- 
ations are too important to simply reflect intrinsic limitations of 
our magnetic model (Sect. 13.2b . they suggest that the large-scale 
field structures seen in polarimetry are more sensitive to dynamo 
fluctuations than the smaller-scale structures that contribute to 
the chromospheric flux. The difference in behaviour between 
the activity tracers and the large-scale field is also visible on 
the shorter term, e.g., between the close-by epochs of 2010.48 
and 2010.59, when a fast evolution of the global magnetic topol- 
ogy (similar in its short timescale to previous observations by 
IPetit et alj|2005h has no clear counterpart in other activity prox- 
ies. 

Several possibilities can be proposed to account for the low 
agreement between the Stokes V and Stokes I activity indica- 
tors. Firstly, Zeeman broadening and chromospheric emission 
are expected to trace a wider range of spatial scales than the po- 
larized Zeeman signatures, which are missing a significant frac- 
tion of the total magnetic flux of the star. Secondly, the differ- 
ent limb visibility of the various proxies is expected to result 
in different temporal evolutions. Thirdly, the flux in chromo- 
spheric indicators tends to saturate with increasi ng magnetic flux 
(ISchriiver et alJ Il989fc lLoukitcheva etal] l2009h . so that strong 



12 



A. Morgenthaler et al.: Magnetic field monitoring of £ Bootis A 



magnetic field variations may result in more gentle Ca II vari- 
ations. 

The complex lon g-term magnet i c varia tions of £, Bootis A, 
previously noted by Ba liunas et al.1 (1 19951) . are reminiscent of 
the behaviour of h i gh-activity, very rapidly-rotating dwarf s 
dDonati et alj l2003bt LTeffers et al] l201lt iMarsden et all 1201 ll) . 
There is obvious magnetic field evolution over the four years 
of observations, but it does not take the simple form of fast 
global pol arity switches observed on other rapidly-rotating sun- 
like s tars dFares et al]l2009fclPetit et alf2 009: Morgenthale r et alj 
1201 ll) . Instead, the main variability observed is a drop of the av- 
erage field strength from 2008 to early 2010, encompassed by 
epochs associated with stronger field. This evolution is also ob- 
served in most of the activity proxies. We note that the toroidal 
field component is relatively weaker whenever the global field 
strength is close to its minimum value, and observe a more ax- 
isymmetric field structure in high-activity states. The magnetic 
topologies of 2007.59 and 2011.07 are very similar, suggesting 
that the repeated observation of this magnetic pattern may be re- 
lated to a type of magnetic cycle that does not imply a global 
polarity switch. Based on our data sets, the large-scale magnetic 
variability of £ Bootis A may seem chaotic, however, we stress 
that not all relevant temporal scales are explored here, so that a 
longer monitoring or a denser temporal sampling may help re- 
veal other patterns. 



5.3. Differential rotation 

Some of the short-term changes in the magnetic topology can 
be taken into account in our model, assuming that they obey 
a simple differential rotation law. The measurement of the ro- 
tational shear was not conclusive at every epoch, presumably 
because of the fast emergence/disappearance of magnetic spots, 
which affects the tracking accuracy of the magnetic tracers used 
to determine the sh ear (a source of error already mentioned 
bv iPetit et alj 12002). Whenever measurable, the values of the 
differential rotation parameters are consistent with a strong sur- 
face shear, with dQ. from five to ten times the solar value. As 
mentioned in Sec. 13.21 part of this observed scatter is likely 
generated by fast changes in the surface topology, but this ran- 
dom effect is not expected to result in systematic biases of dQ. 
measurements. The dQ values are far above some measure- 
ments obtained for very rapidly rotatin g dwarfs in this tem- 
perature domain (see iBarries et al. 2005, who report a drop of 
differential rotation for stars cooler than the Sun). It is, how- 
ever, reminiscent of the strong shear observed on th e T Tauri star 
v2247 Oph, in spite of a spectral type as late as M2 (Donat iet alJ 
120101) . An increase in dQ with rotation is also consistent with 
other observational studies dDonahue et alii 996; Saar 2009]) and 
with the models of lBrown et al.1 (120081) and lBallot et alJ (120071) . 

The high latitudinal shear of £ Bootis A also suggests that 
in spite of the high activity level of the star, Maxwell stresses 
are still inefficient at forcin g solid-body rotation as proposed , 
e.g., for active M dwarfs dMorin et al J 120081: iBrowningl 12 008). 
but a l so for very rapidly rota ting Sun-like stars ( Henry et alj 
119961: iReiners & Sc hmitt 2003]). If such an extreme state is not 
reached by £ Bootis A, the sharp apparent fluctuations of the 
shear level, if real, may indicate that a feedback of the magnetic 
field onto the large-scale flows is probably active in the convec- 
tive layers and is affecting the differential rotation, as suggested 
for oth er high-activity dwarfs (iDonati. Collier Cameron. & Petit! 
l2003al) and from dynamo models of iBrun et al.1 (l2005lh 
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Table 3. Journal of observations for 2007.59, 2008.09, and 
2009.46. From left to right, we list the year of observation, the 
Julian date, the error-bar in Stokes V LSD profiles, and the phase 
of the rotational cycle at which the observati on was made, takin g 
the same rotation period and phase origin as lPetit et al.l d2005l) . 



Year 


Julian date 


(T LSD 


rot. phase 




(2,450,000+) 


10~ 5 / c 




2007.59 


4308.35 


5.4742 


0.8726 




4310.35 


1.7360 


0.1835 




4311.35 


2.1059 


0.3387 




4312.34 


2.7435 


0.4928 




4313.34 


3.8929 


0.6485 




4315.34 


2.8654 


0.9586 




4316.34 


2.7608 


0.1152 




4322.34 


3.2082 


0.0473 




4323.34 


3.7479 


0.2030 


2008.09 


4484.72 


5.9851 


0.3017 




4484.77 


6.3988 


0.3092 




4485.77 


5.2155 


0.4648 




4488.76 


3.2560 


0.9291 




4489.74 


3.4732 


0.0827 




4489.75 


3.5767 


0.0840 




4491.70 


4.0728 


0.3865 




4492.75 


3.3416 


0.5499 




4493.76 


2.7275 


0.7075 




4495.77 


3.4875 


0.0206 




4499.77 


12.7455 


0.6427 




4501.77 


3.5052 


0.9538 




4503.76 


3.1664 


0.2626 




4506.77 


2.9219 


0.7304 




4507.76 


1.3428 


0.8843 




4508.76 


4.5386 


0.0400 




4509.73 


7.8197 


0.1914 




4510.77 


4.1865 


0.3524 




4512.76 


3.1274 


0.6627 


2009 46 


4980 38 


4.0385 


386"! 




4984.45 


12.1589 


0.0192 




4985.43 


3.7381 


0.1722 




4994.47 


5.3269 


0.5786 




4995.46 


4.3773 


0.7318 




5001.44 


14.4977 


0.6617 




5002.36 


10.5647 


0.8056 




5006.56 


3.2475 


0.4579 




5010.46 


5.4779 


0.0654 




5011.44 


4.6241 


0.2175 




5013.55 


4.3254 


0.5450 




5017.38 


6.9613 


0.1415 




5018.40 


3.1676 


0.2992 



Table 4. Same as0for 2010.04, 2010.48, 2010.59, and 201 1.07. 



Year 


Julian date 


O" LSD 


rot. phase 




(2,450,000+) 


IQ- S I C 




2010.04 


5180.78 


8.7507 


0.5527 




5181.78 


3.2145 


0.7082 




5202.71 


4.7321 


0.9641 




5215.65 


4.5171 


0.9763 




5222.71 


3.1795 


0.0743 




5224.69 


3.2308 


0.3823 




5240.62 


4.6358 


0.8593 




5241.66 


7.1186 


0.0211 




5242.63 


3.7747 


0.1718 


2010.48 


5354.42 


3.9811 


0.5580 




5370.45 


4.8856 


0.0514 




5379.41 


3.8653 


0.4445 




5382.37 


3.0234 


0.9054 




5383.37 


2.9261 


0.0600 




5384.37 


2.6600 


0.2154 




5388.36 


2.9740 


0.8374 




5390.37 


5.0528 


0.1487 




5391.40 


2.9108 


0.3093 




5392.40 


2.3126 


0.4650 


2010.59 


5403.41 


3.0789 


0.1766 




5412.36 


3.4796 


0.5697 




5414.35 


2.8979 


0.8791 




5415.35 


2.9610 


0.0335 




5416.36 


4.4882 


0.1919 




5417.35 


3.6313 


0.3448 




5418.38 


3.6560 


0.5047 




5421.36 


3.2580 


0.9690 




5427.34 


2.9531 


0.8985 


2011.07 


5578.71 


3.6620 


0.4409 




5584.76 


3.3521 


0.3814 




5586.75 


3.1925 


0.6901 




5587.75 


3.1746 


0.8459 




5588.76 


2.9784 


0.0028 




5593.75 


4.0095 


0.7787 




5596.71 


4.1288 


0.2387 
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Table 5. Activity tracers of £ Bootis A derived from the Stokes I profiles for each observation for 2007.59, 2008.09, and 2009.46. 



Frac. year 


Julian date 


B, 


Rad. vel. 


Velocity spans 


Line widths 


NcallH 








(G) 


(km.s" 1 ) 


(km.s- 1 ) 


(km.s" 1 ) 






2007.59 


2454315.34 


22.0 ± 0.8 


1.86 


0.02 ± 0.02 


15.059 ± 0.041 


0.4425 ± 0.0005 


0.3577 ± 0.0002 




2454316.34 


8.5+0.9 


1.86 


0.02 ± 0.02 


15.135 +0.038 


0.442 ± 0.0005 


0.3584 ± 0.0002 




2454322.34 


11.7 + 0.9 


1.84 


0.03 ± 0.02 


15.133 +0.041 


0.4427 ± 0.0005 


0.3614 + 0.0002 




2454323.34 


-0.6+1.1 


1.87 


0.03 ± 0.03 


15.385 ± 0.037 


0.4635 ± 0.0006 


0.3631 ±0.0002 




2454308.35 


18.2 ± 1.7 


1.88 


0.06 ± 0.02 


15.444 ± 0.007 


0.4401 ±0.0011 


0.3571 ±0.0004 




2454310.35 


7.3+0.5 


1.88 


0.03 ± 0.02 


14.96 ± 0.032 


0.4419 + 0.0003 


0.358 ± 0.0001 




2454311.35 


1.4 + 0.6 


1.85 


0.04 ± 0.02 


15.352 ± 0.037 


0.4454 ± 0.0003 


0.3585 ± 0.0001 




2454312.34 


2.8 + 0.8 


1.82 


0.07 ± 0.02 


14.875 ± 0.032 


0.438 ± 0.0004 


0.3582 ± 0.0002 




2454313.34 


6.9 + 1.1 


1.84 


0.04 ± 0.02 


14.477 ± 0.037 


0.435 + 0.0006 


0.3566 ± 0.0002 



2008.09 


2454499.77 


7.4 ±3.6 


1.67 


0.05 


± 


0.02 


14.799 ± 0.023 


0.4124 ±0.0018 


0.3542 ± 0.0006 




2454501.77 


9.8 ± 1.0 


1.8 


0.07 


± 


0.03 


14.616 ±0.041 


0.4138 ±0.0006 


0.3562 ± 0.0002 




2454503.76 


4.5 ± 0.9 


1.71 


0.08 


± 


0.02 


14.876 ± 0.034 


0.4309 ± 0.0005 


0.3557 ± 0.0002 




2454506.77 


3.1 ±0.9 


1.73 


0.02 


± 


0.02 


14.137 ±0.036 


0.4039 ± 0.0005 


0.3535 ± 0.0002 




2454507.76 


5.9 ±0.4 


1.75 


0.07 


± 


0.02 


14.467 ± 0.036 


0.4206 ± 0.0002 


0.354 ± 0.0001 




2454508.76 


0.8 ± 1.3 


1.72 


0.03 


± 


0.02 


14.783 ± 0.037 


0.4322 ± 0.0007 


0.3583 ± 0.0003 




2454509.73 


6.7 ± 2.0 


1.73 


0.09 


± 


0.02 


14.97 ± 0.022 


0.4187 ±0.0012 


0.357 ± 0.0005 




2454510.77 


4.6 ± 1.1 


1.68 


0.05 


± 


0.02 


14.533 ±0.036 


0.427 ± 0.0007 


0.3549 ± 0.0003 




2454512.76 


7.1 ±0.8 


1.65 


0.04 


± 


0.02 


13.956 ± 0.024 


0.4145 ±0.0005 


0.3526 ± 0.0002 




2454484.72 


3.9 ±2.1 


1.73 


0.06 


± 


0.02 


14.885 ±0.038 


0.4242 ± 0.0012 


0.3565 ± 0.0004 




2454484.77 


6.4 ±2.1 


1.73 


0.05 


± 


0.02 


14.379 ±0.031 


0.4224 ±0.0012 


0.3573 ± 0.0004 




2454485.77 


3.9 ± 1.5 


1.72 


0.04 


± 


0.02 


14.412 ± 0.038 


0.4212 ±0.0009 


0.3541 ±0.0003 




2454488.76 


1.8 ±0.9 


1.84 


0.06 


± 


0.03 


14.417 ± 0.037 


0.4149 ± 0.0005 


0.353 ± 0.0002 




2454489.75 


5.2 ± 1.0 


1.82 


0.05 


± 


0.02 


15.025 ± 0.038 


0.4221 ± 0.0005 


0.3563 ± 0.0002 




2454489.74 


4.7 ± 1.0 


1.82 


0.04 


± 


0.02 


14.958 ± 0.032 


0.4188 ±0.0005 


0.3562 ± 0.0002 




2454491.7 


-0.4 ± 1.2 


1.75 


0.03 


± 


0.02 


14.38 ± 0.036 


0.4273 ± 0.0007 


0.3555 ± 0.0003 




2454492.75 


-2.3 ± 1.0 


1.73 


0.01 


± 


0.02 


14.495 ± 0.037 


0.411 ±0.0005 


0.3538 ± 0.0002 




2454493.76 


3.7 ±0.8 


1.76 


0.03 


± 


0.02 


14.426 ± 0.042 


0.4079 ± 0.0004 


0.354 ± 0.0002 




2454495.77 


10.5 ± 1.0 


1.78 


0.04 


± 


0.02 


14.619 ±0.037 


0.4314 ±0.0006 


0.3562 ± 0.0002 


2009.46 


2454984.45 


1.9 ±3.5 


1.81 


0.03 


± 


0.02 


14.324 ± 0.029 


0.4098 ± 0.0022 


0.3537 ± 0.0007 




2454985.43 


3.3 ± 1.1 


1.82 


0.04 


± 


0.02 


14.228 ± 0.023 


0.4205 ± 0.0006 


0.3544 ± 0.0002 




2455017.38 


10.4 ± 1.8 


1.88 


0.03 


± 


0.02 


14.423 ± 0.038 


0.4201 ±0.0011 


0.355 ± 0.0004 




2455018.4 


7.5 ±0.9 


1.89 


0.03 


± 


0.03 


14.206 ± 0.018 


0.4252 ± 0.0005 


0.3526 ± 0.0002 




2454994.47 


16.7 ± 1.6 


1.9 


0.06 


± 


0.03 


14.177 ±0.017 


0.4168 ± 0.001 


0.3515 ±0.0003 




2454995.46 


16.1 ± 1.4 


1.9 


0.05 


± 


0.02 


14.467 ± 0.022 


0.4132 ±0.0008 


0.3537 ± 0.0003 




2455001.44 


22.1 ±4.0 


1.89 


0.05 


± 


0.02 


14.017 ± 0.026 


0.4235 ± 0.0036 


0.3521 ± 0.0008 




2455002.36 


9.9 ±2.9 


1.88 


0.03 


± 


0.02 


15.047 ± 0.036 


0.4235 ±0.0015 


0.355 ± 0.0005 




2455006.56 


7.7 ±0.9 


1.89 


0.01 


± 


0.02 


14.409 ± 0.036 


0.4212 ±0.0006 


0.3525 ± 0.0002 




2455010.46 


4.8 ± 1.4 


1.91 


0.03 


± 


0.02 


14.249 ± 0.036 


0.4263 ± 0.0009 


0.3541 ± 0.0003 




2455011.44 


3.3 ± 1.3 


1.9 


0.03 


± 


0.02 


14.271 ±0.037 


0.421 ± 0.0008 


0.3529 ± 0.0003 




2454980.38 


6.5 ± 1.2 


1.85 


0.04 


± 


0.02 


14.287 ± 0.033 


0.4236 ± 0.0007 


0.3534 ± 0.0003 




2455013.55 


-1.8 ± 1.3 


1.88 


0.03 


± 


0.02 


14.265 ± 0.037 


0.4247 ± 0.0009 


0.3512 ±0.0003 
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Table 6. Same as|5]for 2010.04, 2010.48, 2010.59, and 201 1.07. 



Frac. year Julian date 



B, 
(G) 



Rad. vel. 

(km.s- 1 ) 



Velocity spans 
(km.s- 1 ) 



Line widths 
(km.s- 1 ) 



Nr., 



2010.04 



2455202.71 
2455240.62 
2455241.66 
2455180.78 
2455242.63 
2455181.78 
2455215.65 
2455222.71 
2455224.69 



-1.1 ± 1.4 

13.7 ± 1.4 
10.3+ 1.8 
7.3 ± 2.4 
-2.5 ± 1.1 
-0.1 ±0.8 
3.7 + 1.3 
5.6 ± 0.9 
-0.4 ± 1.0 



77 
78 
74 
78 
79 
68 
83 
.8 
84 



0.02 ± 0.02 
0.02 ± 0.02 
0.02 ± 0.02 
0.05 ± 0.03 
0.02 ± 0.02 
0.04 ± 0.02 
0.04 ± 0.03 
0.02 ± 0.02 
0.03 + 0.02 



14.016: 

13.607 : 
13.633 : 
14.057 : 
13.913 : 
13.87 i 
13.779 : 
13.87 i 
13.403 : 



: 0.038 
: 0.034 
: 0.033 
: 0.037 

: 0.017 

0.037 
: 0.038 
0.037 
: 0.018 



0.4093 = 
0.393 ± 
0.3965 = 
0.3983 = 
0.4056 = 
0.4091 = 
0.3997 
0.3921 = 
0.4001 = 



: 0.0008 
0.0008 
: 0.0013 
: 0.0019 
: 0.0007 
: 0.0005 
t 0.001 
: 0.0005 
: 0.0006 



0.353 = 
0.3508 
0.3515 
0.3501 
0.3524 
0.3523 
0.351 = 
0.3502 
0.35 + 



: 0.0003 
t 0.0003 
t 0.0004 
t 0.0005 
t 0.0002 
t 0.0002 
: 0.0003 
t 0.0002 
0.0002 



2010.48 



2455379.41 
2455412.36 
2455382.37 
2455414.35 
2455383.37 
2455415.35 
2455384.37 
2455354.42 
2455416.36 
2455417.35 



12.2: 
10.5: 
10.7: 

4.3 i 

11.2: 

5.2 4 



: 1.1 
:0.8 
:0.9 

0.8 

: 1.1 

0.8 



10.5 ± 1.6 
3.3+0.9 
3.5 + 0.7 
14.4+1.4 



93 
94 
93 
89 
95 
98 
95 
92 



0.07 : 
0.04: 
0.01 : 
0.04: 
0.03: 
0.07: 
0.03: 
0.02: 
0.04: 
0.04: 



0.02 
0.02 
0.02 
0.02 
0.02 
0.03 
0.03 
0.02 
0.02 
0.02 



13.941 ± 
14.318 4 
14.093 : 
14.331 + 
14.128 + 
14.024 : 
13.962 ± 
14.292 ± 
14.023 + 
14.365 + 



0.038 
0.037 

:0.04 

0.026 
0.008 

:0.04 

0.036 
0.038 
0.037 
0.041 



0.4003 : 
0.4115 : 
0.3997 : 
0.3988 : 
0.403 + 
0.4157 : 
0.4025 
0.3967 : 
0.4059 r 
0.398 + 



: 0.0007 
: 0.0005 
: 0.0005 
: 0.0004 
0.0006 
: 0.0005 
t 0.001 
: 0.0005 
: 0.0004 
0.0009 



0.3489 : 
0.3523 : 

0.3511 : 

0.3492 : 
0.3514: 
0.3498 : 
0.3515: 
0.3514: 
0.3512: 
0.3502 : 



: 0.0002 
: 0.0002 
: 0.0002 
: 0.0002 
: 0.0002 
: 0.0002 
: 0.0002 
: 0.0002 
: 0.0003 
: 0.0002 



2010.59 



2455418.38 
2455388.36 
2455421.36 
2455390.37 
2455391.4 
2455392.4 
2455427.34 
2455370.45 
2455403.41 



6.9+1.0 

12.3 + 0.8 

17.4 + 0.8 
17.6+ 1.3 
7.8 + 1.0 
4.3 ± 1.0 
12.8+ 1.0 
10.8 ± 0.9 
10.7 ± 0.9 



92 
91 
97 
97 
95 
92 
92 
93 
95 



0.03 ± 0.02 
0.02 ± 0.02 
0.05 ± 0.02 
0.05 ± 0.03 
0.05 ± 0.02 
0.06 ± 0.02 
0.05 ± 0.02 
0.04 ± 0.02 
0.08 ± 0.02 



14.095 
13.566 : 
13.82 + 
14.429 : 
14.143 = 
14.256 -_ 
13.641 = 
13.723 z 
13.992 = 



t0.04 
: 0.037 
0.039 
: 0.038 
: 0.037 
: 0.025 
: 0.039 
: 0.031 
: 0.033 



0.4111 : 

0.3873 : 
0.3921 : 
0.4058 : 
0.4175 : 
0.4207 : 
0.3889 : 
0.4018 : 
0.399 4 



: 0.0006 
: 0.0004 
: 0.0005 
: 0.0008 
: 0.0006 
: 0.0006 
: 0.0006 
: 0.0005 
0.0005 



0.3522 : 
0.3516: 
0.3496 : 
0.3508 : 
0.3514: 
0.3502 : 
0.3469 : 
0.3537 : 
0.3514: 



: 0.0002 
: 0.0002 
: 0.0002 
: 0.0003 
: 0.0002 
: 0.0002 
: 0.0002 
: 0.0003 
: 0.0002 



2011.07 



2455596.71 
2455578.71 
2455584.76 
2455586.75 
2455587.75 
2455588.76 
2455593.75 



7.3 d 
8.0 4 
6.9 4 
6.9 4 
6.8 4 

15.9: 

7.3 4 



1.1 
0.9 
0.9 
0.9 
0.8 

:0.8 

1.1 



84 
85 
81 
86 
85 



1.8 
1.82 



0.02 ± 0.02 
0.07 ± 0.03 
0.05 ± 0.02 
0.03 ± 0.02 
0.03 ± 0.02 
0.04 ± 0.02 
0.05 ± 0.03 



15.061 : 
14.777 : 
14.489 : 
14.489 : 
14.198 : 
14.667 : 
13.985 : 



0.037 
0.027 
0.025 
0.025 
0.033 
0.032 
0.035 



0.4467 : 
0.4258 : 
0.4325 : 
0.4326 : 
0.4172 : 
0.4311 : 
0.4109 : 



: 0.0007 
: 0.0006 
: 0.0005 
: 0.0005 
: 0.0005 
: 0.0005 
: 0.0007 



0.3615 : 
0.3584 : 
0.3585 : 
0.3541 : 
0.3558 : 
0.3586 : 
0.3542 : 



: 0.0003 
: 0.0002 
: 0.0002 
: 0.0002 
: 0.0002 
: 0.0002 
: 0.0002 



